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Abstract — We present a general problem formulation for 
optimal parameter estimation based on quantized observations, 
with application to antenna array communication and process- 
ing (channel estimation, time-of-arrival (TOA) and direction-of- 
arrival (DOA) estimation). The work is of interest in the case 
when low resolution A/D-converters (ADCs) have to be used to 
enable higher sampling rate and to simplify the hardware. An 
Expectation-Maximization (EM) based algorithm is proposed for 
solving this problem in a general setting. Besides, we derive the 
Cramer-Rao Bound (CRB) and discuss the effects of quantization 
and the optimal choice of the ADC characteristic. Numerical and 
analytical analysis reveals that reliable estimation may still be 
possible even when the quantization is very coarse. 

Index Terms: Quantization, MIMO channel estimation, 
TOA/DOA estimation, EM algorithm, Cramer-Rao Bound, 
stochastic resonance. 

I. Introduction 

In multiple-input multiple-output (MIMO) communication 
systems, where low power and low cost are key requirements, 
, it is desirable to reduce the ADC resolution in order to save 
power and chip area (TJ. In fact, in high speed systems the 
sampling/conversion power may reach values in the order 
of the processing power. Therefore, coarse analog-to-digital 
converters (ADCs) may be a cost-effective solution for such 
applications, especially when the array size becomes very large 
or when the sampling rate becomes very high (in the GHz 
range) [2|. Naturally, this generates a need for developing new 
detection and estimation algorithms operating on quantized 
data. 

An early work on the subject of estimating unknown pa- 
rameters based on quantized can be found in Q. In |4), Q, 
the authors studied channel estimation based on single bit 
quantizer (comparator). In this work, a more general setting 
for parameter estimation based on quantized observations will 
be studied, which covers many processing tasks, e.g. channel 
estimation, synchronization, delay estimation, Direction Of 
Arrival (DOA) estimation, etc. An Expectation Maximization 
(EM) based algorithm is proposed to solve the Maximum 
a Posteriori Probability (MAP) estimation problem. Besides, 
the Cramer-Rao Bound (CRB) has been derived to analyze 
the estimation performance and its behavior with respect to 
the signal-to-noise ratio (SNR). The presented results treat 
both cases: pilot aided and non-pilot aided estimation. We 
extensively deal with the extreme case of single bit quan- 
tized (comparator) which simplifies the sampling hardware 



considerably. We also focus on MIMO channel estimation and 
delay estimation as application area of the presented approach. 
Among others, a 2x2 channel estimation using 1-bit ADC is 
considered, which shows that reliable estimation may still be 
possible even when the quantization is very coarse. In order to 
ease the theoretical derivations, we restrict ourselves to real- 
valued systems. However, the results can be easily extended 
and applied to complex valued-channels as we will do in 
Section rVTl 

Our paper is organized as follows. Section HI1 describes the 
general system model. In Section[IIIl the EM-algorithm operat- 
ing on quantized data is derive and the estimation performance 
limit based on the Cramer-Rao Bound (CRB) is analyzed. 
In Section IIV1 we deal with the single-input single-output 
(SISO) channel estimation problem as a first application, then 
we generalize the analysis to the multiple-antennas (MIMO) 
case in Section [V] Finally we handle the problem of signal 
quantization in the context of Global Navigation Satellite 
Systems (GNSS) in Section [VT] 

Notation: Vectors and mattices are denoted by lower and 
upper case italic bold letters. The operators (») T , (») H , 
tr(»), (•)*, Re(») and Im(») stand for transpose, Hermitian 
ttanspose, trace of a matrix, complex conjugate, real and 
imaginary parts of a complex number, respectively. Im denote 
the (M x M) identity matrix. Xi is the i-th column of a 
matrix X and Xij denotes the (ith, jth) element of it. The 
operator E s i g [»] stands for expectation with respect to the 
random variable s given q. The functions p(s,q) and p(s\q) 
symbolize the joint distribution and the conditional distribution 
of s and q, respectively. Unless otherwise noted, all integrals 
are taken from — oo to +oo. Finally, 1 = xt symbolizes that the 
equality holds for the single bit case. 

II. System Model 

As mentioned before, we start from a general signal model, 
described by: 



r = Q{y), with 

y = f(x,9) + r], 



(i) 

(2) 



where y is the unquantized receive vector of dimension N, 
/(•,•) is a general multidimensional system function of the 
unknown parameter vector 9, to be estimated, and the known 
or unknown data vector x, while tj is an i.i.d. Gaussian noise 
with variance at in each dimension. We assume that the noise 



variance is known, although this part of the work can 
be easily extended to the case where is part of 9. The 
operator Q(-) represents the quantization process, where each 
component is mapped to a quantized value from a finite set 
of code words as follows 

Ti = Q( yi ), if y t e HVH. 0) 

Thereby r\° and r" p are the lower value and the upper limits 
associated to quantized value r$. Additionally, we denote the 
prior distribution of the parameter vector by pe(9) when 
available. Similarly the prior p x (x) is also known, and can 
for instance be obtained from the extrinsic information of 
the decoder output. The joint probability density function 
involving all random system variables reads consequently as 



p(r,y,x,0) = I D{r) (y)- 



1 



(2tt)*< 

where I denotes the indicator function taking one if 



p x (x)pe(8), 
(4) 



y e D(r) = {ye R N \rf < y t < r t up ; Vi € {1, . . . , AT}} , 

(5) 

and otherwise. Note that this special factorization of the 
joint density function is crucial for solving and analyzing the 
estimation problem. A factor graph representation of the joint 
probability density is given in Fig.[T]to illustrate this property. 
Each random variable is represented by a circle, referred to 
as variable node, and each factor of the global function (|4|l 
corresponds to a square, called functional node or factor node. 




(v) 



to- 



Fig. 1 

Factor graph representation. 



III. Construction of the Estimation Algorithm 
and Performance Bound 

Given the quantized observation r, and the log-likelihood 
function 

L{9) = In / / p(r, y, x, 9)dxdy = \np(r, 9), (6) 



our goal is to find the MAP estimate 9 given by 

9 = argmax£(0). 
e 



(7) 



Naturally, the MAP solution 9 has to satisfy 

V e £(0) = 0. 

This condition can be written as: 
V e p(r,0) 



(8) 



p(r,9) 

Vep(r,y,x,0] 



dxdy 



p(r,9) 

Vep(r,y,x,0) p(x,y\r,8) 



p(r,8) p(x,y\r,9 
\7 g p(r,y,x,9] 



dxdy 



p(x,y\r,9)dxdy 

p(r,y,x,8) 
= E !B , ylrt e[Vekip(r,y,x,8)}±0. (9) 

There is also another way to write the optimality condition, 
by first integrating out the variable y to obtain the conditional 
probability of the quantized received vector 



p(r\x,8) = 



1 



||v-/(e,«ou2 



N 



dy 



rr-fi(x,9). 



*(- 



fi(x,8) 



)), 
(10) 

where $(a;) represents the cumulative Gaussian distribution 
reading as 

*(x) = -±=[ exp(-i 2 /2)dt. (11) 

V 27T J-oo 

Therefore, we also obtain an alternative condition as 
V e C{9) = [ VofP^y^dy^ 



p(r,9 

\7 p(r,x,9) ^ p(x\r,9) 
p{r,8) p(x\r,9) 
Vep{r,x,8) 



dx 



p(r, x, 9) 



p(x\r, 9)dx 



V x \ rfi [V e \np(r,x,8)]=0, 



(12) 



which can be explicitly written as, 

2j E a>|r,0 



(e ^ -e ^ )V fi(x,8) 



1-bit 



/ j x\r,i 



VePe(8) 
Pe(8) 

VeP S {8) 
Po(8) 



e 2 < V fj(x,0) 



27rcr T? <&( 



= 



0. 



(13) 



where the last step holds for the single bit case, i.e. r*j € {±1}- 



A. EM-Based MAP Solution 

In general, solving ( PT3T l is intractable, thus we resort to the 
popular Expectation Maximization (EM) algorithm as iterative 
procedure for solving the condition (0 in the following 
recursive way 



,y\r,o l 



Ve lnp(r, y,x,9 



0. 



(14) 



Thus, at each iteration I the following two steps are performed: 

E-step: Compute the expectation 

g(r, 6, 6 ) = ^ x y]r ^ [lnp(r, y, x, 9)} + const 

= E_..,_ * [2y T /(0, s) + ||/(fl, x)||*j/(2^)+lnp e (0) 



E 



x\r,6 



2(ji(o ,x) + E[ v \ x ,r,e }) T f(e, x y 



f(0,x)\\i /(2a 2 rl )+hxp 8 (e), 



(15) 



where 



(r^-J'^.e'))^ 



(>i -/i(«>.8 )) 2 
5^ 



E[^|a!,r,d] = - 



1-bit CTri 

= r, ; 



' — 5^2 — 



/27T ) -v 

M-step: Solve the maximization 

9 = argmaxg(r, 9, 9 ). 



(16) 



In many cases, this maximization is much easier than (0, as 
we will see in the examples considered later. 

B. Standard Cramer-Rao Bound (CRB) 

The standard is the lower bound on the estimation 

error for any unbiased estimator, that can be obtained from 
the Fisher information matrix J (9) under certain conditions 

V{(9-9)(9-9) T } t (JW)- 1 . (17) 
Hereby, the Fisher information matrix reads as [6| 
J=E rl g[V e C(e)Vj£(9)} 



= E. 



r\9 



E x \ r .e [V e \np(r,x,9)} 
Ex|r,e [Vj\np(r,x,9)] 



E n 



\r.6 



(r i L ' t '-/. i (a ! .a))^ (.r-?-fi(°>,Q)V 

^ -e 



E, 



\ r ,e 



(r""-/ 4 (<.,e) 



('■i »)) 2 



Q( r : P -J-M )-$( r .-/;^ 9 ) ) 



v£/i(M) 



(18) 



1 The standard CRB, in contrast to the Bayesian CRB, holds for a 
deterministic parameter, i.e. the prior pg (0) is not taken into account. 



In the pilot-based estimation case (x is known), it simplifies 



to 



(e 



2 V a ft( a! ,e)Vj/ i ( a; ,g) 



1-bit 



E 



2vrcr2$( 



/i(g»,e)' 



d>(- 



(19) 

Additionally, in the low SNR regime (<j ? , > |/i(as,0)|), (fT9l l 
can be approximated by 

o2 



(20) 



where the factor 



1 Y-(e 

PQ = 7T~ 



2tt^ $(r^)_$(i!l) 



< 1, 



(21) 



depends only on the quantizer characteristic (here assumed to 
be the same for all dimensions) and represents the information 
loss compared to the unquantized case at low SNR, in the pilot- 
based estimation case. For the single bit case, i.e., (r l °, r up ) G 
{(— oo, 0), (0, oo)}, the Fisher information loss pq is equal to 
2/tt, which coincides to the result found in [7|, |8| in terms 
of the Shannon's mutual information of the channel. For the 
case that we use a uniform symmetric mid-riser type quantizer 
[9 1, the quantized receive alphabet is given by 

ri G{(-|--i + fe)A; fc=l,... ,2 b } = K, (22) 

where A is the quantizer step-size and b is the number 
of quantizer bits. Hereby the lower and upper quantization 
thresholds are 

r lo fn-f for ri >-#(2>-2) 
] — oo otherwise, 



and 



up 



A 

2 



for n < - 
otherwise. 



-(2 b -2) 



In order to optimize the Fisher information at low SNR ( f20b 
and get close to the full precision estimation performance, we 
need to maximize pq from OTT) with respect to the quantizer 
characteristic. Table U shows the optimal (non-unform) step 
size A opt (normalized by a 2 ) of the uniform quantizer de- 
scribed above, which maximizes pq for b £ {1,2,3,4}. If 
we do not restrict the characteristic to be uniform, then we 
get the optimal quantization thresholds which maximize pq 
in Table [TTJ We note that the obtained uniform/non-uniform 
quantizer optimized in terms of the estimation performance is 
not equivalent to the optimal quantizer, which we would get 
when minimizing the distortion, for a Gaussian input J9). l n 
addition, contrary to the quantization for minimum distortion 
the performance gap between the uniform and non-uniform 
quantization in our case is quite insignificant, as we can see 
from both tables. 



Table I 

Optimal Uniform Quantizer. 



b 


A op t 


PQ 


1 




2/vr 


2 


0.704877 


0.825763 


3 


0.484899 


0.945807 


4 


0.294778 


0.984735 



Table II 

Optimal non-Uniform Quantizer. 



b 


Optimal thresholds 


PQ 


1 





2/tt 


2 


0;±0.704877 


0.825763 


3 


0;±0.306654;±0.895 107 ;± 1 .626803 


0.956911 


4 


0;±0.143878;±0.4204768;±0.708440;±1.017896; 
± 1 .364802;± 1 .780058;±2.346884 


0.989318 



In the following, the theoretical findings will be applied to 
the channel estimation problem and for a GNSS problem with 
quantized observations. 

IV. Example 1: SISO channel estimation 

We first review the simple problem of SISO channel esti- 
mation considered in 0. 

A. Pilot-based single bit estimation (one-tap channel) 
The SISO one-tap channel model is given by 

r t = sign(y 4 ) = sign(/iXi + m), for i G {1, . . . , N}, (23) 

where N is the pilot length and xi € { — 1, 1} is the transmitted 
pilot sequence with normalized power. The channel coefficient 
h € R is here our unknow parameter, i.e. 6 = [h]. 
It can be shown by solving the optimality condition ([TJt 
(with uniform prior pe(6)) that the ML-estimate of the scalar 
channel from the single bit outputs r^, is given by 



h 



(24) 



Besides, the Fisher information (T% becomes in this case 
J(h) 



Ne 



27Rj2$( 



)*(- 



(25) 



This expression of the Fisher information is shown in Fig. [2] 
for N = 200 as function of h 2 /a 2 . In Fig. |3] the CRB, i.e. 1/ J 
and the relative exact mean square error (MSE) of the ML- 
estimate from TV = 200 observations, both normalized by h 2 , 
are depicted as function of the SNR= h 2 /a 2 . We interestingly 
observe that above a certain SNR, the estimation performance 
degrades, which means that noise may be favorable at a certain 
level, contrary to the unquantized channel. This phenomenon 
is known as stochastic resonance, which occurs when dealing 
with such nonlinearities. We can naturally seek the optimal 
SNR that maximizes the normalized Fisher information, i.e. 
minimizes CRB/h 2 : 

h 2 

= argmax- ; ,"Z, ; , — = 2.4807 = 3.9458dB. 



O p fc 



27T 7 $(77)$(- x /7) 



This results obtained by numerical optimization of the Fisher 
information coincide with the results found in [5| through 
observations at asymptotically large N. 




4„ *5 6 

h I at (linear) 



Fig. 2 



Fisher Information vs. cri for a SISO channel, b = 1, N = 200. 




2 3 4 

h 2 /cr^ (linear) 



Fig. 3 



Estimation error vs. <A for a SISO channel, 6 = 1, N = 200. 



(26) 



B. Pilot Based Estimation (two-tap channel) 

Now let us consider a more general setting with a two-tap 
inter-symbol-interference (ISI) channel 

n = sign(?/i) = signfioXi+hiXi^+rn), for i G {1, . . . ,7V}, 

(27) 

where ho and hi are the channel taps. Again we utilize a 
binary amplitude pilot sequence x E { — 1,1}^ and we try to 
find the the ML-estimate of the parameter vector 8 = [ho, hi] T 
in closed form. Ignoring the first output r\, the ML-condition 



( fT3] l turns to be (pg(6) = 1) 



C. Non-Pilot Aided (Blind) Estimation 



N 



= 0, 



N 

E 

i=2 



1 " 



(h Q x i + h 1 x i _ 1 ) 



ri(hoXi-\-hiXi-i) ^ 



(28) 



Suppose now that a unknown binary symbol sequence Xi £ 
{+1, .1} is transmitted over an additive white Gaussian noise 
(AWGN) channel with an unknown real gain h. The analog 
channel output is 



y t = h- Xi + r], 



(33) 



Taking the sum and the difference of these equations delivers 
respectively 



N 



C^O^i + ^l^i — l) 

72- 



Ee i 
rAxi + Xi-i) 77 -t >— = 0, 



t=2 
JV 



(29) 



r<(a;i — ) 77 rr r- = 0. 



i=2 



Next, we multiply the numerator and denominator of each 
equation by g(- r '( h °*«+ h *"-i? ) to get 



2_^n(xi + x l - 1 ) ( ho+hl ) , fe n+ /T7 T 

i=2 V u a „ > 



i=2 
N 



(30) 



= 0. 



' S)( jjo—hl. \S)( ha — h-i_ \ 

Then, using the fact that 

2$(-r^i + hl ) = 1 - rix 4 erf(^-i— 
where erf(-) denotes the Gaussian error function, we get 

N N h 

/Xxj + Xi-x)ri = (N -1 + 22 XiXi-i)eTt{—i= ] 

i=2 i=2 V2(Tr) 



(31) 



i=2 
JV 



i=2 
JV 



^2(xi - Xi-i)ri = (N - 1 - XiXi-i)eTf( h ° 1 ). 

i=2 i=2 V 2 ' 7 ') 

(32) 

Finally, solving the last equations with respect to ho and hi, 
we get the ML solution. 



( 



ff- 1 



■ N 

i=2 



erf" 1 
V 



AT 

N+X) a;»a;j_i 

i=2 

■ TV 

i=2 



rf" 1 



• N 

J2(xi-Xi-i)r t 

i=2 



iV 



+err 1 



JV 
i=2 

V N 

J2(xi-i-Xi)ri 

i=2 



N 



N+J2 XiXi-x 



) 
) 



The solution consists in quite simple computations (apart of 
the final application of erf -1 ) since we only have to do with 
binary data {r^Xi € {±1})- 



where the variance of the noise rj, a^, is also unknown. Ad- 
ditionally, the receiver is unaware of the transmitted symbols 
Xi. Based on N quantized observations i\ = Q(yi), we wish 
to estimate the parameter vector 6 — [h, cr r) ] T . We note an 
inherent ambiguity in the problem: the sign of the gain h and 
the sign of Xi cannot be determined individually. We also note 
that at least 2 bits are needed in this case, because a single 
bit output does not contain any information about h. Since 
the ML problem is intractable in closed form, we resort to the 
EM approach. The EM-update for h can be obtained from the 
general expressions in ([T5T > and ( TTo*l > as 



N ^ 



:,a:e{+l,-l} 



fc«[*(!i^)-*(S^L)]-^L(e 



.1.2 



*(^)-*(^) + <T>l 



r]°-h l • 



-.1,2 



TT7T 



— e 



+ c 



" (34) 

while the update for the noise variance follows from the 
expectation 



n 



V2al 



N 



h l )e 



2 v \ - 

lo_rl,2 



•r? p -W - 



(35) 



where we used the conditional distribution 



p{r)i\r u a l h 1 )-. 



y2ir<j 7 ; Y27TCT,, 



A°-hK 



(^6) 



The Cramer-Rao Bound that can be easily obtained from the 
likelihood function 

C(h,a v ) = 



In $(■ 



-*c- 



Jo . 



as well as the MSE of the estimates h and a. 



(37) 
found by 



Monte Carlo simulation, both normalized by h , is depicted 
in Fig. |4] as function of h 2 /a 2 , where the pilot length is 

N 

the MSE of the estimate a v also exhibits the non-monotonic 
behavior mentioned before with respect to the SNR. 



100 and the quantizer resolution is b = 3. Clearly 



CQ 
« 

id 

CO 

5 





-e-MSE-h ' 




MSE- o 

1 

-e-CRB-h 




CRB- a 

'1 










-~ e --oT^ 





5 10 15 20 25 30 35 

h 2 /a^ (linear) 
Fig. 4 

MSE AND CRB OF THE BLIND ESTIMATES OF h AND t7 v VS. h 2 /<T% FOR A 
SISO CHANNEL, b = 3, N = 1000. 



V. Example 2: Pilot-based MIMO channel 
Estimation 

Now, we consider the MIMO case. We begin first with the 
problem of estimating a 2 x 2 MIMO channel from single bit 
outputs, since it can be also solved in a closed form, as shown 
later on. 

A. Single bit estimation of a 2 x 2 MIMO Channel 

As example, let us consider a pilot-based estimation of a real 
valued 2x2 channel matrix assuming a single-bit quantizer 

Ti = sign(/iii£Ci + h l2 x 2 + 7]i), (38) 

where x\, x 2 E {—1,1}^ are the pilot vectors transmitted at 
each Tx antenna, while r\, r 2 E {— 1, 1}^ are the received 
vectors at each Rx antenna. The maximum likelihood (ML) 
channel estimate 6 — [hn, hi2, hi, hix] T can be found by 
solving (fT3l in closed form, similarly to the 2-tap SISO 
channel (see Subsection HV-Bt , as 




(aji + x 2 ) T r t 



N + x{x 2 



-erf 



Xj) T r t 



N - x{x 2 



(39) 



with i, j E {1,2}, j = 3 — j and erf -1 the inverse function 
of the error function. We can see that the HW implementation 
of the estimation task is still considerably simple, since only 
shift registers, counters and a look-up table for erf -1 would 
be necessary. Fig. [5] shows the Monte Carlo simulation and the 
CRB of the estimation error J^- E[(/ijj — hij) 2 ] for a given 
2x2 channel as function of the noise variance a^. Thereby, 
the Fisher information matrix can be obtained from ( fl9] > as 

x];x 2 )e 



J = 



47TCr2$(^a)$(- 



hi+h 2 ' 



(N - xjx 2 )e 



(h 1 -h 2 y 
— 



(40) 



As example, we took the specific channel matrix 



H 



2 

0.5 



1.5 
-1 



(41) 



Fig. [5] shows also the MSE when using the analog (unquan- 
tized) output, which is exactly 



MSEI 



[Xi 



x 2 ] T [x 1: 



*2])" 



(42) 



Clearly, the estimation error under quantization does not 
increase monotonically with higher a 2 , as we know already 
from the SISO case. 



0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




Unquantized 

- - 1-bit CRB 
1-bit simul. MSE 



10 



n 

Fig. 5 

Estimation error vs . a* for a 2 x 2 real valued channel, 
N = 200. 



B. Pilot-based MIMO channel estimation of arbitrary size 

Let us consider now a more general setting of a quantized 
linear MIMO system 

y = vec[HX'} + rj, (43) 

and 

r = Q(y), (44) 



pMxiV 



with a channel matrix H e M LxM and X' e 
contains N pilot- vectors of dimension M. Hereby, we stack 
the unquantized, quantized and the noise signals into the 
vectors y, r and rj, respectively. Our unknown parameter 
vector is therefore 8 = h = vcc[H] and we have the system 
function 



f(h,X) = Xh, 



(45) 



oM-NxM-L 



where the new matrix X £ I J "'" AJ "'" contains again 
the pilot-vectors in a proper way. Furthermore we assume, 
contrary to the previous cases, that a priori distribution p(h) 
is given according to h ~ Af(0, Rh)- With this definition the 
EM-iteration <JT3J and ( fTol l reads in this case as 



E-step: Compute for i = 1, . . . , N 



27T 



$( ^-w« )_$( tM*y]i ) 



(46) 




-5 ,5 

-101og 10 (^) 

Fie. 6 



Estimation error vs. <r^ for a real valued 4x4 MIMO channel, 
6 = l, N = 1000, R h = lie, £i, 3 6 {-1, +1}. 



M-step: 



h l+1 = (X T X + a 2 ^ 1 )- 1 ■ X T (Xh! + b l ). (47) 

Let us at this point validate the convergence of the EM- 
algorithm to a unique optimum solution. For this, we write 
the log-likelihood function explicitly 



£(0)=5>(*( 



(48) 



This log-likelihood function is a smooth convex function with 
respect to 9. This follows from the log-concavity of 



<&(& - z) — $(a - z), 



(49) 



b > a, with respect to z, since it is obtained from the convolu- 
tion of the Gaussian density and a normalized boxcar function 
localized between a an b, which are both log-concave [10|. 
Therefore, the stationary point of the EM-iteration fulfilling 
the condition ( fT3l is the unique optimal solution. 
Fig. [6] illustrate the average MSE defined by 



MSE = E 



h-h 



(50) 



under different bit resolution for a 4x4 MIMO channel with 
i.i.d. unit variance entries. Hereby, we chose an orthogonal 
pilot sequence, i.e. X T X — Rh — Iw with Xij G { — 1, +1}- 
The estimation error in the unquantized case, which is given 
by 



MSE 



6— >c 



<#r (( 



X T X 



(51) 



is also shown for comparison. Obviously, at medium and 
low SNR, the coarse quantized does not affect the estimation 
performance considerably. 



VI. Example 3: Quantization of GNSS Signals 

The quality of the data provided by a GNSS receiver 
depends largely on the synchronization error with the signal 
transmitted by the GNSS satellite (navigation signal), that is, 
on the accuracy in the propagation time-delay estimation of the 
direct signal (line-of-sight signal, LOSS). In the following we 
will study the effect of quantization in terms of simulations and 
the CRB as derived in Section IIII-BI We assumed an optimal 
uniform quantizer as given in Table [TJ We will first assess the 
accuracy of a standard one-antenna GNSS receiver in case no 
multipath is present. Secondly, we will assess the behavior of 
array synchronization signal processing in a multipath scenario 
applying the innovative derivation of the EM algorithm as 
shown in Section IIII-AI This assessment is based on the work 
presented in IfTTII . In the following we assume a GPS C/A code 
signal with a chip duration T c = 977.52 ns, a code length of 1 
ms and a bandwidth of B = 1.023 MHz. The received signal 
is sampled with the sampling frequency f s = 2B. We only 
use one code period as an observation time where the channel 
is assumed constant during this observation time. 

The synchronization of a navigation signal is usually per- 
formed by a Delay Lock Loop (DLL), which in case no mul- 
tipath signals are present, efficiently implements a maximum 
likelihood estimator (MLE) for the time-delay of the LOSS t\. 
In Fig. [7] the lower bound of the RMSE of t\, for different the 
number of bits, is given in terms of the J CRB(fi) in meters 
where Co denotes the speed of light. A nominal SNR for a 
GPS C/A signal is approximately —20 dB. In Fig. [7] one can 
observe that the yj CRB(fi) does not significantly decrease 
further for more than 3 bits, thus a rather simple hardware 
implementation is sufficient for such a GNSS receiver. 

Now, we assess the EM algorithm as derived in Section 
IIII-AI with pe(0) being a uniform distribution, hence consid- 
ering a ML estimator. We consider a two path scenario where 
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Fig. 7 

RMSE OF f ■ CO VS . BIT RESOLUTION AND SNR WITH ONE ANTENNA 

f s = 2.046MHz. One code period is used for estimation. 



the LOSS and one reflective multipath signal are received by 
an uniform linear antenna array (ULA) with M = 8 isotropic 
sensor elements. We define 



e = [M7l T ,M7}VV T i T ] T 



(52) 



12 3 4 

b 

Fig. 8 

RMSE OF fl ■ C AND f 2 ■ C (IN METER) VS. BIT RESOLUTION FOR 
M = 8, 0i = -30°, 02 = 62°, At = 0.3T C , SNR = -22.8DB, 
SMR = 5dB, Au = 0Hz. One code period is used for estimation. 



with the vector of complex amplitudes 7 = [7i,72] T , the 
vector of azimuth angles <j> = [0i,(/>2] T , the vector of time- 
delays t = [ti,T2] t , and the vector of Doppler frequencies 
v = [vi,V2] T . The parameters with the subscript 1 refer 
to the LOSS and parameters with the subscript 2 refer to 
the reflection. The reflected multipath and the LOSS are 
considered to be in-phase, which means arg(7i) = arg(72), 
and the signal-to-multipath ratio (SMR) is 5dB. Signal-to- 
noise ratio (SNR) denotes the LOSS-to-noise ratio and we 
assume SNR = -22.8dB. The DOAs for the LOSS and the 
multipath are <f>\ = —30° and fa — 62° respectively. Further, 
we define the relative time-delay between the LOSS and the 
multipath as At = |n — T2I = 0.3T C and relative Doppler 
Av = \vi - v 2 \ = 0Hz. In Fig. E the RMSE of fi and f 2 
vs. the bit resolution is depicted. In Fig. [9] the RMSE of <fii 
and 4>2 vs. the bit resolution is shown. Based on the results 
presented in Fig. [8] and Fig. [9] one can derive the important 
statement that 4 bits seem to be sufficient for high-resolution 
estimates with respect to the considered channel conditions. 

VII. Conclusion 

A general EM-based approach for optimal parameter esti- 
mation based on quantized channel outputs has been presented. 
It has been applied in channel estimation and for GNSS. 
Besides, the performance limit given by the Cramer-Rao 
Bound (CRB) has been discussed as well as the effects of 
quantization and the optimal choice of the ADC characteristic. 
It turns out that the gap to the ideal (infinite precision) case in 
terms of estimation performance is relatively small especially 
at low SNR. This holds independently of whether the quantizer 
is uniform or not. Additionally, we observed that the additive 
noise might, at certain level, be favorable when operating 



3.5 



1—. 3 








m 




K 


2.5- 







^> 




w 


2- 






1 




s 


1.5- 



R MSEfji) 
RMSE(<£ 2 ) 



0.5 

1 Z 3 4 

b 

Fig. 9 

RMSE OF 01 AND 02 (IN DEGREE) VS. BIT RESOLUTION FOR M = 8, 
02 = -30°, 02 = 62°, At = 0.3T C , SNR = -22.8DB, SMR = 5dB, 
Au = 0Hz. One code period is used for estimation. 



on quantized data, since the MSE curves that we obtained 
were not monotonic with the SNR. This is an interesting 
phenomenon that could be investigated in future works. 
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